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FOREWORD 


This  report  supersedes  Air  Materiel  Conunand  Memorandum 
Report  MCREXA5-4595-8-2,  "A  Method  for  Predicting  Dynamic 
Landing  Loads",  prepared  by  Lee  S.  Wasseman  under  date  of 
20  February  194s.  The  purpose  of  re-writing  this  report  is 
to  expand  and  revise  details,  to  present  an  additional 
derivation  of  the  theory  and  to  present  a new  comp\itation 
form. 


This  report  was  prepared  in  the  Dynamic  Loads  Section, 
Dynamics  Branch,  Aircraft  Laboratory,  Directorate  of 
Laboratories,  Wright  Air  Development  Center  under  Research 
and  Developnent  Project  136?,  Structural  Design  Criteria, 
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ABSTRACT 


This  report  supersedes  Memorandum  Report  MCHEXA5-459 5-8-2, 

”A  Method  for  Predicting  Dynamic  Landing  Loads”,  20  February 
1948.  The  purpose  of  rewriting  is  to  correct  minor  erix)rs  and  make 
several  refinements. 

Dynamic  responses  may  be  computed  as  the  sum  of  the  rigid 
bocfy  response  and  the  vibratory  responses  in  each  normal  mode. 

The  rigid  bocy  re^onse  is  determined  first  from  basic  airplane 
parameters  and  in  this  report  is  assumed  trapezoidal  in  ^ape. 

This  traoezold  is  then  applied  to  Ihe  equation  of  motion  of  the 
elastic  ayst«n  to  determine  the  vibratory  response. 

The  vibratory  response  of  an  elastic  system  to  a trapezoidal 
forcing  function  can  be  computed  algebraically  or  graphically.  The 
algebraic  computation  method  is  motivated  by  two  distinct  principles; 
discontinuity  and  stperposition;  and  the  graphical  computation 
method  is  motivated  by  the  superposition  principle.  New  computation 
forma  are  provided  for  both  the  algebraic  and  the  graphical  methods. 


Three  particular  problems  ar«  solved  to  compare  theoretical 
and  measured  results,  to  serve  as  a computation  guide  and  to 
illustrate  the  flexibility  of  the  ^proach.  In  the  first  problem 
the  effect  of  varying  basic  parameters  is  discussed  with  a flow 
chart. 
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INTRODUCTION 


The  eTolublon  of  rastly  stepped  xq)  perfomance  of  airplanes 
has  placed  Increased  Importance  on  deslgnixig  structures  to  close 
tolerances  to  minimize  mei^t  penalties.  The  problem  of  pre- 
dicting and  analyzing  landing  loads  provides  a fertile 

field  for  replacing  ea^irieal  criteria  by  rational  ooBqputatlons. 

This  report  presents  an  acceptable  method  for  the  conputatlon 
of  dynamic  landing  loads.  The  basic  assumption  of  a trapezoidal 
shape  for  the  rigid  bo^y  loads  was  suggested  by  Mr.  Lee  Wasserman. 
This  assumption  is  in  agreement  with  test  results.  Fortunately  it 
is  also  simple  to  handle  as  a forcing  function  for  the  differential 
equation  of  motion  of  the  vibratory  system. 

The  present  paper  supersedes  MCHBXA5-4595->d-2,  saae  title » 
dated  20  Februaiy  1948*  The  object  of  this  revision  is  to: 

1.  Correct  the  following  errors  in  the  examples: 

Page  in 


Reads 

Should  Read 

13,  26  ami  33, 
column  10 

n 

O)  n 

33,  column  12, 
row  "A"  to  "B" 

oo 

— oo 

These  errors  had  little  effect  on  the  final  results. 

2.  Expand  and  revise  the  details  of  the  original  theory 
for  easier  reading. 

3*  Provide  an  additional  derivation  of  the  acceleration 
response  to  a trapezoidal  forcing  function* 

4*  Present  a new  self-contained  computation  form  >bieh 
eliminates  a significant  amount  of  superfluous  arithmetic. 

5.  Present  a flow-chart  analysis  of  the  basic  parameters 
in  the  first  example.  The  effect  of  changing  basie  parameters  is 
shown  to  be  intricate  but  predictable. 
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In  the  examples  of  Sections  III  through  V experimental  data 
were  used  to  determine  fuselage  and  landing  gear  frequencies  and 
the  relerant  modes  of  vibration.  For  other  aircraft  es^erlmental 
values  of  frequencies  may  not  be  available.  The  computation  of 
fuselage  frequencies  presents  no  problem  (see  Section  IV);  but 
the  computation  of  landing  gear  frequencies  requires  further  investi- 
gation. 


If  It  Is  not  known  which  natural  vibration  modes  are  relevant, 
calculatlms  may  be  necessary  beginning  with  the  mode  of  lowest 
frequency  and  continuing  until  the  cooputed  responses  are  no  longer 
significant.  Reference  9 discusses  the  difficult  question  of  select- 
ing relevant  modes. 
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SECTION  I 


BASIC  THEORY 


1.  It  is  assumed  that  at  the  moment  of  landing  wing  lift 
exactly  counterbalances  the  weight  of  the  aiiplane,  so  that  the 
vertical  velocity  is  constant.  Consequently  the  response  is  due 
entirely  to  dissipation  of  kinetic  energy  at  impact  and  may  be 
considered  in  two  stages: 

a.  Rigid  bo<fy  response  of  the  whole  structure. 

b.  Vibratory  response  within  the  structure  in  each 
normal  mode, 

2.  Energy  equilibritim  conditions  must  be  satisfied: 

a.  For  the  structure  as  a whole  (dete mining  rigid  body 
response)  i.e,.  Kinetic  Energy  »=  Potential  Energy. 

b.  For  each  particle  of  mass  (detemining  vibratory 
response)  i.e..  Inertial  Work  - Elastic  Work  = External  Work, 

3.  a.  The  rigid  body  response  is  dete  mined  from  the 
equilibrium  conditions  for  the  structure  as  a whole.  Briefly, 
the  kinetic  energy  at  impact  is  determined  from  the  rate  of 
descent  and  gross  weight.  But  this  is  equal  to  the  potential  energy 
of  the  tire  and  strut  work.  The  tire  deflection  vs  load  curves 
detemine  the  tire  work,  and  the  strut  work  is  detemined  assuming 
isothermal  expansion  and  quasi-adiabatic  compression. 

Assuming  a trapezoidal  shape  for  rigid  body  load,  the  time 
history  is  now  easy  to  evaluate.  For  further  detail  of  this  method 
see  Reference  3, 

b.  The  vibratory  response  in  eadi  mode  is  detemined  by 
the  local  equilibrium  conditions: 


X Inertia  work  - work  = work  (1) 

In  particular,  let  the  entire  mass  be  represented  by  a finite  number 
of  elements  W/  (e.g.  gear,  boc^,  tail,  wing,  etc.)  each  located  at 
a point  in  space.  Then  the  vertical  displacement  and  acceleration  of 
the  element  ttt/  in  each  mode  dqoend  on  the  element's  location  and  can 

be  written  C/X  and  c/x  respectively,  where  Cf  is  a constant 
determined  by  the  position  and  mode. 
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Now  consider  incremental  displacements  the 

elements /»/•  caused  by  external  forces  where  is  the 

time  history  of  the  external  forcing  function  with  unit  amplitude. 
Ten5)orarily  neglecting  damping: 


Force 

X 

Distance 

* Work 

Inertia;  m/CfX 

X 

Ci  efx 

« fm,  Ci  x)(ci  <Jx) 

• 

Ela  Stic : A/  Ci  X 

X 

a <Jx 

» (Hi  a xUci  tJx) 

External;  ffF(t) 

X 

Ci  Jx 

. {/i  frtXci  <Jx) 

And  so  equation  (1)  can  be  written: 

^ Frt)  (2-a) 


or 

(lmiCf^)k^(XKic))x  ^(Zfjc,)‘f(f)  (2-b) 


It  will  now  be  shown  that 

particular  case  (Zfffi  Cf  *Jx  * CZKi  cp )x 

and  there  is  sinple  harmonic  motion  so  that:  Ji^-cu^X 
is  the  mode  frequency.  Thus  (Z  X7jCi*')  ^Z  Ki  c-^ 

Substituting 

(ZmiCi*‘)x  ^(Zmici*’)ai^x*  (ZfiCi)'/^M 


For  in  the 
where  ui 


(2-c) 


or 


**  z 
X ^ CU  X 


/T(f) 

ZmiCj*- 


(2-d) 


4*  The  effect  of  structural  danping  can  be  approximated  using 
a dimensionless  "damping  coefficient"  y acting  on  the  displacement 
so  that: 

X ^ CU f y ; X * ^ 

Z mi  a 

It  is  this  equation  (3)  which  gives  the  vibratocy  acceleration  response 
X . In  this  report  ^ is  rather  arbitrarily  assianed  to  be  .10, 

5.  It  is  sometimes  convenient  to  by-pass  the  (constant) 
coefficient  of  defining: 

Genarallzed  force 

GAF  > Generalized  Acceleration  factor  « Generalized  mass  (4) 

^In  "g"  units.^  - 

______________________  a ■* 

* /T/  “ ^ring  constant 
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6,  This  eqmtion  (3)  can  be  solved  separately  for  eadi  noraal 
mode  because  of  the  original  assumption  that  nomal  modes  do  not 
feed  energjr  to  each  other  or  to  the  rigid  body  modes.  The  assumption 
is  validated  by  the  excitation  of  reasonably  pure  natural  modes  during 
ground  vibration  tests, 

7,  The  next  problem  is  to  determine  which  normal  modes  are 
important,  (See  Reference  9 for  further  discussion,)  The  experience 
of  AMC  dynamic  tests  shows  that  only  the  first  few  modes  are  important 
unless  there  is  appreciable  coupling  between  landing  gear  fore  and 
aft  vibrations  and  hi^er  structural  modes, 

8,  An  outline  of  the  computational  procedure : 

Step  I;  Ccanpute  rigid  body  vertical  load  time  history  from 
basic  airplane  parameters. 

Step  II ! Compute  rigid  body  drag  load  time  history  if 
appropriate.  Assume  (empirically)  the  coefficient  of  friction  is  ,55 
until  the  wheel  gets  to  speed,  after  which  the  rigid  bo(ty  drag  load 
falls  to  zero  in  one-quarter  the  spinup  time. 

Step  III;  Determine  which  modes  of  vibration  are  important. 
Experimental  data  of  frequencies  will  be  used  if  available;  If  data 
is  not  available  equations  of  motion  of  the  structure  can  be  used, 
b\jt  care  is  required  in  selecting  the  appropriate  degrees  of  freedom. 

Step  IV;  Compute  generalized  acceleration  factor  if  appropriate. 

Step  V;  Compute  vibratory  response  in  each  mode.  The 
theory  for  this  computation  is  discussed  in  Section  II, 

Step  VI;  Obtain  the  time  history  of  total  acceleration 
or  structural  force  by  appropriate  combination  of  rigid  and  vibratory 
components.  Notice  that  the  trapezoid  considered  as  rigid  body 
component  may  have  a different  ordinate  from  the  trapezoid  considered 
as  forcing  fimction  for  vibrations. 
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SECTION  II 

DERIVATION  OF  THE  ACCEUERATION  RESPONSE  TO  A 
TRAPEZOIDAL  FORCING  FUNCTION 
(Solution  of  Equation  (3)) 

1.  Given  forcing  function  * 


2,  It  is  required  to  determine  explicitly  the  vibratory 
acceleration  response  for  each  interval  AB,  BC**..In  symbols,  find 
y e:q}licltly  where  x for  each  interval  is  given  implicitly  by: 


with  T defined  as  follows: 

AB:  T*t  and 

7$ 

(5) 

BC: 

7i 

CD: 

0^  7"^  n 

(6) 

DE: 

EF: 

0*  rs 

— F: 

0£  /"*" 

* The  function  has  beim  sketched  in  sufficiently  general  terns  to 
satisfy  all  common  problems.  In  some  caaes,  e»g.,  and 
(pos.  or  neg.). 
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fora: 


3.  It  is  apparent  * that  (5)  has  an  explicit  solution  of  the 

• .7-.  y -ri  (7) 


X » & 


foAmcuT-t-Jt  coAoiTJ 


For  example  in  AB  the  initial  conditions  at  7*^  are  X‘0  and 
Then:  X = -4:^  (3) 


oi 


The  prolaLem  is  thus  solved  in  the  interval  AB. 


4.  To  get  the  solution  in  each  successive  interveil  two 
approaches  may  be  taken: 

a.  Discontinuity  Derivation:  Evaluate  the  constants 
and  Jk  in  (7^  by  applying  the  initial  conditions, 

b.  Stiperposition  Derivation:  Resolve  trapezoid  into  the 
sum  of  straight  lines  all  of  which  emanate  from  points  on  the  t-axis. 

Each  of  these  approaches  tdll  now  be  demonstrated.  Of 
coiirse  they  must  lead  to  identical  results. 


Discontinuity  Method 
(Paragraphs  5-9) 

5.  It  is  required  to  determine  the  constants  a and  ^ in  (7) 

in  each  interval  from  the  initial  conditions  Xo  and  , (The 

values  of  x and  x when  T'-tf.) 

6,  The  solution  in  AB,  according  to  paragraph  3#  is  (8),  Thus 
the  value  of  x and  x will  be^known  at^  the  end  of  AB.  Similarily 
once  the  interval  BC  is  solved,  x and  x will  be  known  at  the  end 
of  BC.  So  if  each  interval  is  solved  in  turn  it  is  fair  to  assume  in 
all  cases  that  X and  x are  known  at  the  end  of  the  previous  interval. 


* For  since  x ui* (/^^J)x^/f7')  is  to  be  solved  for  x (not  x ), 
in  present  form  it  is  an  integral  eqmtion  (involving  x which  is 
).  To  get  a differential  equation  in  x differentiate  (5) 
twice  with  respect  to  7*  to  get 


X ^ X •o 

since  ///V  is  always  linear  in  T so  that  /(T)*0 

This  is  a standard  hrao^neous  equation  in  X with  ths  general 
solution:  (<»f»f>rox}maf0) 


X - 


la  sfn  Oi 


T 


cos  Ui 
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7.  ^ According  to  paragraphs  5 ^d  6,  it  will  be  sufficient  to 
evolve  Xo  and  Xo  from  and  Xm  (where  'e*  denotes  evalxiation 
at  the  end  of  the  previous  interval).  Since  x is  continuous  from 
one  interval  to  the  next, 

Xc-X^  (9) 

But  X is  not  continuous,  so  that  Xo  ^ Xe. 

So  a gimmick  will  be  introduced  to  evaluate  Xo  • Basically, 
we  transport  the  discontinuous  variable  x through  the  medium  of 
the  continuous  variable  X T 

First  consider  equation  (5)  when  differentiated  once: 

(10) 

Evaluating  (30  at 

Xtt  = ^ •"  ) . 

From  continxiity,  Xo  « X0 

So  evaluating  (10)  at  the  end  of  the  previous 

***’*  f / -^ ) X^  * 

or,  from  (12);  ^ f J ) - jtV 

Finally,  from  (11)  and  (H), 

(15) 

8,  Solving  (7)  and  its  derivative  at  T*  O, 

and  applying  (9)  and  (15),  the  constants  and  can  be  (ietermined; 

N j 

X«  - 

Xo  * *»UJ  - 


a ^ %r 


<0  (11) 

(12) 

interval; 

(13) 

(U) 
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9.  In  particular,  the  solution  in  each  interval  is; 


a.  In  AB: 


A' 


S/ncut 


at  B: 


//z 
~ to  T, 


ai  7i/z 


S/ntuTt 


(18) 


(19) 


Xb 


r. 


^/>7a;  5^(20) 


b.  In  BC: 

Xo^Xa  Z^Xb--^  (21) 

77 

So  from  (17) 

//a  ^ sMte 

too  ZtoTi 

a.  co&co^^  ^ _dk.  (22) 

oj  Ti  ioTJ" 


-9  to  Tt/i.  y-  _ 

<7  » ^ & /cos  tu7^  ~ ^ s/a  u 


to 


r, 


So  that  (7)  becomss: 
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This  may  b«  rewritten: 


X « 


(dj/* 


(24) 


where 


I 


cos  <o 


ojTi 


» af'c  Am 


/P4 


d*Z 


(25) 

(26) 


Note: 


The  choice  of  the  correct  quadrant  for  is  crucial. 

The  quadrant  must  be  selected  so  that: 


has  the  algebraic  sign  of 


i*o*:  %ir\ipf  has  the  algebraic  sign  of 
and  has  the  algebraic  sign  of 


(27) 


Finally  at  C: 

Xc  * sy/itUiTi  *^t) 

Xc.  • & ^ i/x^  .3^  s/nfdf/i  i*^)/'Af«5/Si»^i‘^28) 


c.  In  CD:  X»  • Xc 


49^  999 

Xo  ^ Xc  ^ - 


from  (15) . 


7k 


(29) 


So  from  (17) 

yi  ^ i/ Xs  CCS  (ou  Tjt  ^*)  (30) 

cuTk 

And  (after  simplification): 

X * .£*  //fjf**  ) (31) 
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where 


^ g castcuTz  ¥-^4) 

€u7s 


<kz  » ‘V'C  /i?/7  Xc 
/^z 


(32) 


Note:  ^ must  be  detemlned  as  follows: 

(a)  Quadrant: 

C Xc 

(cas  ^zj  algebraic  sign  of  ^ /p^ 

(b)  Angle: 


(33) 


/i/7?  *^Zm  Xc 


Finally 


Xi»^&  ^WcoT^r^^z) 

d«  In  D£  the  process  is  similar,  resulting  in: 

<$>3) 


(34) 


(35) 


CQACtu^'/'^z)— 

where:  \ ^ 

<^a  » a/c  /cm  X4> 

(36) 

Note:  ^must  be  detemined  as  follows; 

(a)  Quadrant:  r»  ■> 

^^Has  the  algebraic  sign  of  (37) 

(b)  Angle:  Asfn  , X» 

XTs 

•t  /^H  37  . 

X^  = G / Xa  /X^3  &//i  »*  ^3^  (3®) 

e.  And  in  EF: 


and 


H 

X e 


-s<^t 


(39) 
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^9^  /Tr% 

vbmT9  ^ g /^3  C0»f<u7^  ¥■  ^m) 

Of  7a 


I 


*t*9  . arc  fan  Xg 

Note:  yvmust  be  determined  as  follows: 

(a)  Quadrant: 

/ssjn  ^9)  S' X9'  P 

cas^^j  the  algebraic  sign  of  ( /Pa  \ 

H 

(b)  An^e:  fan  ^ , Xs 


(40) 


(a) 


and  A>  » 


-S*** 


7S>^ 


fPa  ^ sfn{uj7k  / 


(42) 


As  stated  in  paragraph  4j  the  derivation  of  paragraphs  5-9 
could  be  replaced  by  paragraphs  10-13 • 


Smperposition  Method 
(Paiagraphs  10-13 ) 

10.  According  to  paragraph  3 > it  is  ea^  to  find  the  response 
to  any  straight  line  emanating  from  the  time  axis*  The  forcing 
function  of  figure  1 can  be  resolved  Into  such  lines: 


Fig.  2.  Linear  Decomposition  of  the 

Trapezoidal  Forcing  Function  of 
Figure  1 
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tdwre  the  slope 


/Ur)  is 

r, 

/Ur) 

77 

/Ut) 

U,'Mz 

n 

/^m 

/7,-7/z 

73 

/Ur) 

_ 

Ts 

/Ur) 

(43) 


Tff 

11,  Then,  according  to  paragraph  3»  the  response  to 


/rr) 

is  S//?ior 

caTi 

t^o 

/zrr) 

•*  , ^9*»»  , t 

Xa  , - tfz  ff  ^ S//?ajr 

a/Tt 

^(r) 

X j . cu  t ‘ 

«/  7k 

/iio 

/^rr) 

/7*-7/z  c ^ s/nuit 

ut  7k 

''  r^Q 

/zrr) 

n 

'•••  ft.0 

12.  But  since  , , ^ ^ r-  ✓ ^ 

/ " / "^/t  */z  -^/z  ^/s 

response  to  /O’)  nust  be  the  sum  of  the  responses  to  / 

(44) 


(45) 


etc. 
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13*  That  la: 


a.  In  ABi 


and 


b.  In  BC:  X*X,-t  Xz 


X - //*  sy^uf/ 

(46) 

tuTf 

Xm»-7ui  e s/noiTt 

(47) 

cu77 

X , //*  ^ sJneo/  (4d) 

*o7/  ctfTt 


But  since 

ft  /“ tj  -jfuiTt/z 

icuf, 


X ~ ^ 


Expanding  and  collecting: 

lUr, 

■^/ Ml  s^ncuTt  ) cos  <ofl 

and  finall7:  ' ^ j 

i.e,  *.) 

0*77  *o7t 

«• 

s arx:  fon  Xa 
A?0 

Note:  ^^iBuet  be  determined  as  foUovs: 

(a)  Quadrant 


fs^n  ) 
I cos  ^7) 


has  the  algebraic  sign  of  y 


(49) 


(50) 


(51) 


(52) 


(53) 
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§9 

(b)  Angle;  /a/^  - Xa 


And 


Xc  a & 


i^Xa  ' s}7>tcu^  »*  ^i) 


(54) 


e.  In  CD 

M #«  M #•«.  •*  **  ** 

X = Xt  i-  Xz  /-  Xa  • fX»  Xa)  Xa 


(55) 


!• 

X 


cuTa 

Bub  since  /'* /”  VTi 

rXa'f'/^f  ^ M-Ajt  js/nm. 

Expanding  and  collecting; 

»*  /^ r~^‘“^'^/^J */!‘' sin/iiTz  *■  *j)  cosmt*7 

od  finallr;  -J 

(59) 


and  finally;  ^ 

;f  = e /xTT^  <f.z) 


i«faere 


— y*^  _ a I ^ . 

” ^ w X^  ^<u7i  ^ Xt  "“//z 

„ (60) 

“ Xcfrf  _Xc_ 

/9i 


Note;  mast  be  determined  as  follows; 
(a)  Quadrant; 


has  the  algebraic  sign  of 


(61) 
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(b)  Angle:  fart^x  » 


d.  In  DB:  J<  *('X.  ^J(x 


/ 

OnTs 

But  since  / "»  /"*  y*  Th 
by  tbe  sane  process  as  before: 

x . /BT^  ^7,rmt  ■ V *3)  (*5) 

wlwre  a /Xc  »*•/?/' <r4»s/2i^^ 

( 6P^  /44N 


/ ^3  » Xx» 

( 

Note:_iA3a(a8t  be  determined  as  foUeers: 

(a)  Quadrant: 

has  the  algebraic  sigi  of 

j CCS 

i4 

(b)  Angle:  /isn  = Xe* 

A^s 


A>s  ^ 

e»  In  W, 


^s<x>  /ji%  ^/f/^  :i//7  ^u/7^  ^ 


X m ( X*  Xz  y*  X3  y-  } y*  xs 


•« 

X » ^ 
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Again  since  /*  “'  * /"  ""  7^ 

_ 

X ^/xS ^ >Pv- ' ■'■> ■**)  (71) 


Where 


^ /xi  COS  Ta/Tif  »* 


cfy^c  /a77 

/^<» 


(72) 


Note:  must  be  determined  as  follows: 


(a)  Quadrant: 


has  the  algebraic  sign  of 


(73) 


(b)  Angle: 


A»r? 


/«"4» 


The  results  cf  paragraph  13  (that  is,  A in  each  interval 
and  at  the  end  of  each  interval)  are  identical  with  the  results  of 
paragraph  9.  From  these  the  zeros,  peaks  and  discontinuities  of 
the  acceleration  response  are  easily  found  and  plotted,  A self- 
explanatory  form  for  the  cumbersome  computations  will  be  used  in 
Sections  III  - V. 
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SEC?ncai  III 


VERTICAL  INCREMENTAL  ACCELERATION  OF  THE 
WING  TIP  OF  AN  F-80A  AIRPLANE  WITH 
FULL  WING  TIP  TANKS 


This  problem  uses  ths  parameters  of  landing  2,  fli^t  37> 
of  the  AMC  F-SOA  tests  reported  In  reference  4« 

Ccmiparisons  (but  not  computations)  are  also  shown  for  the 
eases  of  half-empty  (landing  30>1)  and  empty  (landing  28-2)  tanks 
in  Figure  3* 

These  particular  landings  were  selected  because  their  values 
of  the  rigid  body  incremental  acceleration  correlated  better 
than  any  others  with  the  theoretical 



U. 

Step  I.  Vertical  Load  Time  History:  ^ 

1.  Basic  airplane  data: 

Gross  weight  ■ 14,000  lbs 

Vff  • Rats  of  descent  • 6 ft/ see  (assumed) 

//  > Mass  per  main  gear  ■>  217  slugs  (2  wheel  landing) 

tv  - Static  load  per  main  gear  • 6250  lbs  (3  point  position) 

/(£  ■ Kinetic  energy  per  main  gear  ■ 3906  ft-lbs  (2  vdieel 

landing) 

ou  ■ Natural  frequency  of  wing  « 16,75  rad/ sec 
2»  Basic  oleo  data: 

Sr  * Total  extension  plus  latent  air  columnN^  > (7*95  + 
2.00)  in.  - .8292  ft 


* 

Latent  air  column  - trapped  air  volume 

piston  cross  section  area 
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^5  - Static  txtonslon  plus  lateit  air  eolusa  • 

(3.00  ♦ 2.00)  in.  - .4167  ft. 

Assuaing  Isothennal  expansion  from  static  to  fully  extended 
position,  the  "load  factor"  at  total  extension  is: 


7? . _ (74) 

Then  assuming  quaAi-adiabatie  compression  from  the  fully 
extended  position  daring  impact,  the  extension  at  any  load  factor 

77  is: 


idiere  7>  ■ load/6250  lbs  and  Y » 1.3  (from  reference  8). 


Basle  tire  data:  < 

[Manufacturer* 

s data ) 

© 

@ 

0 

Incremental  Tire  Work  <■ 

Tire  Deflection 

Tire  Work  * ■*'0 

Load  (Lbs) 

(Ft)  _ _ 

(Ft-Lbs) 

(Ft-Lbs) 

2500 

.058 

73 

73 

6500 

.125 

302 

375 

9000 

.166 

318 

693 

Total  work: 

@ @0 
BZ92®  0-0 

@ 

@ 

77^^ 

6ZBO 

.40 

1.356r 

® 

\ 

1.04 

.4832 

.5715 

.4285 

.3553  2309  2684 

1.44 

.3490 

.4455 

.5545 

.4598  4138  4831 

• (average  Talue  of  load  during  increment)  (deflection  in 
increment); Approximating  the  area  \mder  the  tire  curve  trapeeoids. 


The  load  is  not  yet  sufficient  to  ooBQ)ress  the  strut. 
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5.  Tire  deflection,  oleo  deflection  and  load  when  kinetic 


energy  per  gear  - total  woric: 

Total  Work 
(Ft-Lbs) 

Load 

(Lbs) 

Tire  Deflections 
(Ft) 

Oleo  Deflection 
(1^) 

2684 

6500 

*125 

.355 

- 3906 

-7923 

Zr  - .148 

Zo  - .415 

4831 

9000 

*166 

*460 

6*  The  time  for  tire  conpresaion  7r  oleo  compression  To 
and  tire-oleo  ej^nsion  Tar  are  determined  from  the  formulas  of 
reference  3»  assuming  a trapesoidal  time  history  for  the  strut  axial 


load: 

(77) 

1 = a./S/ 

If 

1 3 a 

|f 

(78) 

Step  II*  Drag  Load  Time  History: 

The  drag  load  is  not  used  since  fore  and  aft  forces  do  not 
put  appreciable  energy  into  the  first  tmcoupled  bending  mode* 


Step  III;  Inqjortant  Modes: 

Since  this  example  involves  the  vertical  acceleration  at  the 
elastic  axis,  the  torsional  mode  of  the  wing  is  e^qsected  to  have 
little  effect*  But  if  the  dynamic  torque  were  to  be  predicted  the 
torsional  mode  and  the  torque  caused  by  the  drag  load  would  be  con- 
sidered* 
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step  IV;  Generalized  Acceleration  Factor  for  Vibratory  Re^onse: 


Reference  4 gives  the  computed  frequencies  and  mode  shapes  for 
the  first  bending  mode  for  full  tanks,  half  full  tanks,  and  empty 
tanks,  GAF  (per  "g"  load  at  gear)  ■ (g  units) 

Where : 

/y  B static  load  per  iidieel,  assuming  all  loads  taken  by 
the  main  gear 

, « vertical  deflection  of  wing  at  gear 

" vertical  deflection  of  wing  tip  (first  bending  mode) 

cj/Af/  = mass  of  wing  element  at  station  "i”  (slugs) 

^ ^ veirtical  deflection  of  wing  at  station  ”1” 

' vertical  deflection  of  wing  tip  (first  bending  mode) 

The  results  per  "g”  landing  = •-1.03  (full) 

- -1.06  (half -full) 

« -1.15  (empty) 

Step  7.  Vibratory  Acceleration  Response: 

1,  Table  1-a  computes  the  vibratory  response  by  the  desk- 
calculator  method.  The  form  is  designed  to  make  the  cuntiersome 
computations  as  mechanical  and  well-grovqsed  as  possible.  It  is  meant 
to  be  self-ersplanatory  after  following  the  theory  in  Section  II. 

2,  Table  1-b  compvites  the  vibratory  response  by  the  graphical 
method. 

Step  VI.  Total  Acceleration  Time  History: 

Figure  3 shows  the  total  tip  acceleration  as  the  sum  of  the 
rigid  body  and  vibratory  accelerations.  The  results  for  half-full 
and  empty  wing  tip  tanks  are  also  shown. 

The  measured  and  computed  results  agree  well  for  full  tanks,  but 
for  empty  tanks  additional  modes  should  probably  be  included  in  the 
computations. 


WADC  TR  54-28 


19 


WADCTR54-28 


20 


GRAPHICAL  SOLUTION  OF  VIBRATORY  ACCELERATION  RESPONSE  TO  A TRAPEZOIDAL  FORCING  FUNC 


^ ^ o o o o 

3 B di  k " 

PPPV 


.-j  ij 

.iviviviv>  S 

\ + I 

X O 3 ^ w > 

. ‘SL  w w M 
+ 3 3 3 3 

:X  + + <■  +* 
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step  VII.  Variation  of  Basic  Parameters: 


1,  In  this  example  many  parameters  were  involved  - about  twelve 
either  previously  measured  or  assumed  and  the  others  computed.  Will  it 
be  possible  to  predict  (without  completely  new  computations)  the  effect 
on  the  final  results  caused  by  changing  any  of  these  parameters? 

2.  To  facilitate  this  discussion,  a flow  chart  (Table  1-c)  of 
the  computations  is  shown.  Each  parameter  is  oriented  vertically  by 
its  order  of  convenient  computation  and  horizontally  according  to  the 
ntnber  of  previous  parameters  >diich  are  used  to  compute  this  one. 

Lines  emanating  upward  from  any  block  lead  to  all  of  the  other  parameters 
which  effect  this  one.  The  chart  gives  little  if  any  quantitative 
help  in  particular  cases. 

3 • But  the  folloviing  general  conclusions  seem  worthwhile : 

a.  Change  in  a basic  parameter  has  an  intricate  effect, 
and  in  many  cases  the  best  procedure  is  to  recompute  from  scratch. 

b.  Thete  are  a few  parameters  ^ S , &A^)  which  effect 
only  the  vibratory  response.  The  effect  of  altering  any  of  these 
is  rather  easily  predicted. 

c«  The  problem  is  conveniently  considered  in  three  stages 
(separated  by  broken  horizontal  lines  in  Table  1-c). 

(1)  Equating  energies, 

(2)  The  rigid-bo(ty  trapezoid. 

(3)  Vibratory  response  and  summation, 

4,  If  many  changes  were  contemplated  in  the  first  section,  it 
may  be  profitable  to  solve  the  second  and  third  sections  only  once 
in  general  terms,  then  plug  in  the  particular  numbers  computed  from 
Part  1.  See  reference  10  for  some  computations  on  alteri^  basic 
parameters. 
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TABLE  1-C:  PARAMETER  FLOW-CHART  (EXAMPLE  l) 
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SECriOH  17 


TAIL  BOOM  INCRIIlEin'AL  ACCELERATIONS 
OF  AN  F-61  AIRPUNE 


This  problem  vises  the  parameters  of  landing  eighteen 
of  the  AMC  tests  reported  in  reference  5*  This  particular  land- 
ing is  selected  because  the  maxianm  rigid  body  vertical  load  is 
closest  to  the  computed  value* 

Step  1*  Vertical  Load  Time  History: 

1*  Basic  airplane  data: 

Gross  weight  • 25,000  pounds  (app. ) 

A/  a Mass  per  main  gear  ■ 3dd  slugs  (2  wheel  landing) 

« Rate  of  descent  ■ S ft/see  (assumed) 

JV  * Kinetic  energy  per  gear  12,400  ft-lbs 

2*  Basie  oleo  data: 

• Total  extension  * ■ 10  in*  ■ *S33  ft 

• Static  extension  * ■ 2.87  in.  ■ *239  ft. 

Assuming  isothermal  esqjjanslon  from  static  to  fully  extended  position, 
the  "load  factor"  at  total  extension  « 

/Jr  m ■ s ^.2B7  (79) 

Then  assuming  quasi-adlabatie  cosqpresslon  from  the  fully  extended 
position  during  impact,  the  extension  at  any  load  factor  "n"  during 
impact  is: 

where  /7  ■ loa<V^12500  lbs  and  7 ■1*3* 


* Assuming  latent  air  column  > 0 
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3.  Basic  tire  data: 

© 


@ 


Tire 

Deflection 

Load  (Lbs) 

(Ft)  _ 

6000 

.125 

13000 

.250 

22000 

.333 

Total  work: 

@ 

© @ 

/2sao  nrQi 

@>^•3 

.43  .5979 

.6733 

1.04  .2760 

.3715 

1.76  .1631 

.2473 

@ 

0 

Incremental 

Tire  Work  - 

Tire  Work  * 

^ 0 

(Ft-lbs) 

(Ft-Lbs) 

375 

375 

1188 

1563 

U53 

3016 

0 0 

@ 

© 

' ^ .833m  ® 

rarr^A 

0-0 

.3267  .272 

1644 

2019 

,6285  .523 

6199 

8362 

.7522  .627 

13794 

16810 

5.  Tire  deflection,  oleo  deflection  and  load  when  kinetic  energsr 
per  gear  ■ total  work: 


Total  Work 
(Ft-Lbs) 

Load 

(Lbs) 

Tire  Deflection 
(Ft) 

Oleo  Deflection 
(Ft) 

8362 

13000 

.250 

.523 

TT^  -0.2400  17302 

Zr  - .290 

Zo  - .573 

16810 

22000 

.333 

.627 

6.  The  times  for  tire  compression  To  , oleo  compression  Tr 
and  tire-oleo  escpanslon  Tor  are  determined  from  the  formulas  of 
reference  3* 

7r  ^ 1^_ 

TWofx  T y(  Tiwtx / (81) 


7a  - 


\ 


2/^Xo  ' p S€^C. 

y 


(d2) 

(33) 


♦ (^  ■ (Average  load  during  increment)  (deflection  in  increment) 
A^roximating  the  area  under  the  tire  curve  by  trapedoids 
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step  II.  Drag  Load  Time  History 

1.  Basic  parameters: 

I.  ■ moment  of  inertia  of  eadi  landing  gear  rolling 
assenbly  = 12,3  slxig  ft^ 

Rjj  « rolling  radius  of  wheel  ■ 1,65  ft. 

= landing  speed  = 147  ft/sec. 

Tire  coefficient  of  friction  = 0.55  (assimed) 

2.  Computation: 

^ = angular  velocity  of  wheel  after  (64) 

spin  up  ^ » 147  rad/sec  ■ 89»1  rad/sec 

1.Q5 

angular  velocity  after  tire  conpression  ■ (85) 

a.ss 7r  * 24.0  rad/sec, 

2Jtjr 

Assimiing  (empirically)  that  peak  drag 
load  ■ ,55  Pfliax 

» angular  velocity  required  for  spin  up  during  (86) 
oleo  compression  « O ~ &r  “ 65.1  rad/sec 

7s  ■ duration  of  skid  during  oleo  compression  ■ (87) 

Oo  JTa  «.0510  sec. 

- total  spin  up  time  = ,0886  sec, 

TJr  e time  to  drop  to  zero  drag  load  » (88) 

* ,0221  sec.  (Empirical  formvila)  ^ 

Step  III,  Important  Modes 

1,  The  vibratory  acceleration  of  the  tail  boom  is  caused 
partially  by  the  vertical  load  and  partially  by  the  drag  load.  It 
is  assumsd  here  that  the  vertical  load  acts  directly  to  set  up  a 
vertical  vibration  in  the  tail  at  the  natural  frequency  of  the  tail. 
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2.  In  order  to  detennine  the  important  vibratory  modes  set  up 
by  the  drag  load,  the  aiiplane  is  assumed  to  have  the  configuration 
shewn  in  figures  4 and  5 with  the  three  degrees  of  freedom  illustrated, 

^ . i / 

^ 


S! 


r^yjL  (mt) 


Fig,  4*  Geometric  Representation 


z<-h} 


Fig,  5.  Positive  Direction  of 

Generalized  Coordinates 


Then  (Lagrange's  equation  for  zero  external  torque) 


Si  s ^ js 


With: 

(90) 

^ *^Z 

for  which  I is  the  moment  of  inertia  of  the  body  and  at,  and  Wjt 
are  natural  frequencies  of  tail  and  gear  reqaectively. 

Directly  from  (^)  and  (90) 

t'7 miJi" /■  rrji,  ^ ^ frft  AJ  /yfz  IzZ  ^ ^ 

«• 

//7,y  y /r>,ui*y  m^J,^  * S (91) 

/y?z  B f 777z  COz  z.  ^ /Til  Jz  ^ ‘ O 

These  equations  are  unwieldy  because  they  cont^n  both  the 
parameters^  , , 2 and  their  second  derivatives  » y » JS  • 
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But  the  spring  action  is  assumed  single  harmonic  motion.  Then: 

where  <3>f  and  are  constants  and  au  is  the  response  frequency  of 
the  system.  The  object  is  to  determine  the  value  of a/  for  each  of 
the  possible  modes  of  action. 

Directly  from  (92): 

Jr^-cu*y  (93) 

Rewrit ihg  (91)  by  substituting  (93)j 

7f2,  //7z  J.B  ‘a 

j,  ^ /■//  ‘ ^ (94) 

This  system  of  three  equations  in  the  coordinates  f3  ^ y , 3^ 
has  all  constant  terms  zero.  Therefore  it  has  only  the  trivial 
solution  B - y = ^ ^ 0 unless  the  determinant  of  the  coefficients 
^ ty  i ^ is  zero.  That  is: 

//>  /»,// 

A O I (95) 

A a [/-(^)] 

The  only  parameter  laiknown  in  (95)  is  to  . For  the  F-6l  airplane, 
the  other  parameters  have  the  values  (for  the  complete  airplane) 

ft^  slugs 

= 6^  slugs  » 700  slugs 

3 9 


COt  **  Wi  " 47.1  rad/sec 
^ = 24  ft  A ^ 7.7  ft 


Expanding  (95)  and  solving  for  to  gives  two  solutions,  and 
corresponding  values  of  the  ratios  B , Z . These  are: 

~T  y 

Ist  Coupled  Mode  2nd  Coupled  Mode 

CU  a CA)i  « -^7.  / /'creZ/s^. 

^/y-0 

2/j  . -J»J'  =-z.e^5  Z j, 

mi.  Zt  ~y  ‘ -jr  ^ 
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step  IV.  Generalized  Acceleration  Factors. 

1,  For  the  vertical  load: 

ff/rf  , f-Aiso) -/=e*4.(97-a) 

per  foot  of  tail  deflection. 

2.  For  the  drag  load: 


a.  In  the  first  mode: 

(/r3QZ\ Q > ^ a 

per  foot  of  tail  deflection. 

b.  In  the  second  mode: 


(97-b) 


Tail  deflection  « 1 f t - ^ 

Thus  gear  deflection  = .3208  ft. 


And  ^ ^ 

So  that: 


-.QZk91  rad. 

(O.S6')  (/7^oz\ 


i'aa^oe) 


02^97)^  y- 


f 2^ 

per  foot  of  tail  deflection 


5.a8-Tj 

(97-c) 


Step  V.  Vibratory  Acceleration  Response. 


1.  Since  the  mode  for  the  vertical  load  has  the  same  frequency 
as  the  first  mode  for  drag  load,  the  two  rigid  boc^r  forcing  functions 
can  be  combined  into  one  trapezoidal  forcing  function  as  shewn  in 
Figure  6. 


Forcing  Functions 
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The  re^onse  to  this  function  is  computed  in  table  2-a. 

2,  The  response  in  the  second  mode  due  to  drag  load  is  computed 
in  table  2-b. 


Step  VI.  Total  Acceleration  Time  History. 

1,  The  total  acceleration  is  the  sum  of: 

a.  Rigid  bo(fy  response  to  vertical  load.  This  and  all 
other  components  of  the  final  response  are  nondimensionalized  in 
units  of  1/2  airplane  iieighb.  Thus  the  maximum  value  of  the  rigid 
body  response  will  be  scaled  as  t73o2.  _ / 

/J2  SOO  “ ' • 


b.  Rigid  body  response  to  drag  load.  The  maximum  value 
is  deterained  from  the  definition  of  torque: 

Torque  - (.55)  (17302)  (7.7)  ft-lbs 

But:  Torque  * I»C  > (1.11  x 10^)fMax.  vertical  accelerationy^24]ft  lbs 

2 

So: 


Max  vertical  acceleration 


3.169  g (98) 


c.  Vibratory  responses  to  vertical  and  drag  loads. 

2.  All  the  above  components  are  summed  and  conqpared  >dth 
experimental  results  in  figure  7, 
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TABLE  2-b 


k 


SECTION  V 


B-17G  LANDING  GEAR  DRAG  LOAD  (NORMAL  TO  STRUT) 


Parameters  are  taken  from  landing  number  1 of  flight  3 of 
the  B-17G  landing  load  tests  of  AMG,  reported  in  reference  6, 

This  example  is  a particularly  good  criteria  of  tiie  accuracy 
of  the  predictions  based  on  the  trapezoidal  theory,  because 
was  actually  measured* 

Step  I.  Vertical  Load  Time  History: 

1,  Basic  airplane  data: 

Gross  weight  ■ 48,137  pounds 

//  = Mass  per  main  gear  ■ 750  slugs  (two  wheel  landing) 

* Static  load  per  gear  ■ 21,775  lbs  (3  point  position) 

■ Rate  of  descent  = 7 ft/sec  (measured) 

« Kinetic  energy  per  gear  =18,380  ft-lba 

oi  ■ Natxiral  frequencgr  of  gear  fore  and  aft  * 73.2  rad/sec 
(measured  from  landing  records) 

2.  Basic  oleo  data: 

“ Total  extension  + latent  air  column  • (9.6  + 1.56)  in,  ■> 
.930  ft, 

« Static  extension  -f  latent  air  column  - (1.44  1.56)  in.  ■> 

.250  ft. 


Assuming  isothermal  expansion  from  static  to  fully  loaded  position, 
the  "load  factor"  at  total  extension  » 


^ /fe  . €P,Z6S8 


(99) 


Then  assximing  quasi-adiabatic  compression  from  the  fully  extended 
position  during  impact,  the  extension  at  any  load  factor  "n"  during 
impact  is: 


(100) 
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iAier9/7  ■ load/21,775  lbs  and  t ■ 1*3  (assumed). 
3.  Basie  tire  data: 


0 

® 

Incremental  Tire  Work  ■ 

Tire  Deflection 

Tire  Work  ♦ 

Load  (Lbs)  (Ft) 

(Pt-lbs) 

(Ft-Lbs) 

10000 

.160 

800 

800 

15000 

.223 

792 

1592 

20000 

.283 

1050 

2642 

25000 

1219 

3861 

Total  Work: 

® 

0 0 

0 

0 @ 

>7  _ (Z) 

Zt77/r 

a9zo-(^  ©'© 

.4592 

.5854  .8339 

.1661 

.1545  1545 

2345 

.6889 

.3902  .4849 

.5151 

.4790  7185 

8777 

.9185 

.2927  .3886 

.6114 

.5686  11372 

14014 

1.1481 

.23a  .3273 

.6727 

.6256  15640 

19501 

5«  Tire  deflection,  oleo  deflection  and  load  whan  kinetic  energy 
per  gear  <■  total  work: 


Tota  Work 
(Ft-Lbs) 

Load 

(Lbs) 

Tire  Deflection 
(Ft) 

Oleo  Deflection 
(Ft) 

1404 

20000 

.283 

.569 

*=  18380 

- 23978 

jCr  - -326 

- .614 

19501 

25000 

.337 

.626 

6.  The  times  for  tire  conpression  7r  oleo  coopression  To 
and  tire-oleo  e^anslon  Tor  are  determined  from  the  fonnulas  of 
reference  3« 


rr  . 


2^  y^Xr' 


tc  s^c (loi) 


* ^ s ( average  load  during  increment)  (deflection  in  increment) 
(^^roximating  the  area  under  the  tire  curve  by  trapezoids) 
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^./960  Atre. 


(102) 


n 


s 


-f-Tr) 


c>.  ^ 97a  s^c. 


(1C3) 


Step  II.  Drag  Load  Time  History: 


!•  Basic  parameters: 

X/t  * moment  of  Inertia  of  each  landing  gear  rolling  assembly 
31.23  slug  ft2 

« rolling  radius  of  %dieel  - 1 .942  ft . 

->  landing  speed  a 147  ft/sec. 

Tire  coefficient  of  friction  ■ .55  (asstimed) 


2 . Confutation: 

^ a angular  velocity  of  vdieel  after  (104) 

spin  up  a a 147  rad/sec  • 

75.70  rad/sec. 

^7  a angular  velocity  after  tire  compression  a (105) 
RkAss a 19,68  rad/sec  (peak 

drag  load  a ,55  ?„,.,). 

a angular  velocity  for  spin  if  during  oleo  (IO6) 

confression  a 0 -- ■ 56.CXi  rad/sec 


7%  a duration  of  skid  diudng  oleo  compression  ■ (107) 

- .0683  sec, 

7t  V’Tjs  a total  spin  up  time  a ,1167  sec, 

7z  « time  to  drop  to  zero  drag  load  a Tr  t**  (108) 

= ,0292  sec.  (empirical  formula).  ^ 
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step  III»  Important  Modes: 


The  first  mode  only,  (frequency  ■ natural  frequency  of  gear  fore 
and  aft  ■ 73.2  rad/sec)  is  assumed  to  be  axcited  by  the  vertical  load 
component  normal  to  the  strut  and  by  the  drag  load  component  normal  to 
the  strut. 

Step  IV.  Generalized  Acceleration  Factor  of  Vibratory  Response: 

1.  It  will  first  be  shown  that  it  is  not  necessary  to  compute 
the  vibratory  response  to  the  normal  component  of  vertical  load  and  the 
vibrating  response  to  the  normal  component  of  drag  load  separately. 
Rather  the  entire  vibratory  response  can  be  computed  from  a "double" 
trapezoid  combination  of  the  two  forcing  functions.  Consider: 


Fig.  8.  Time  Histories  of  Vertical 
Load  and  of  Drag  Load  in 
Units  of 


Fig.  9.  Load  Normal  to  Strut  at  Time  t 

Due  to  Vertical  Load  and  Drag  Load 


WADC  TR  54-28 


3a 


2.  Thus  the  normal  to  the  strut  due  to  rigid  body  motion  is; 

^ ^y/?  S • “*  ^ ^ ^ (109) 

3*  It  is  therefore  possible  to  consider  a combined  function  to 
represent  rigid  body  acceleration  and  as  forcing  function  for  the 
vibratory  acceleration: 


Fig.  10.  The  Composition  of  the 

Trapezoidal  Forcing  Function 


4.  In  Exao^les  I and  II  ^e  vibratory  acceleration  X and 

rigid  body  acceleration  •X'/>5>are  found  and  added  to  find  the 

total  acceleration  response. 

5.  But  in  this  exaaple  the  structtural  force  rather  than  the 
total  acceleration  is  to  be  computed.  That  is,  in  this  example  find 
structural  force  normal  to  the  strut  at  the  ajde  due  to  a one -foot 
vertical  deflection  of  the  axle*  Writing  equation  (1)  with«6r  a/ft; 


c*  t^nfi)'x  y-  (' (HO) 

where:  mi) Is  required  structural  force 

(S^ci^mj)  X is  the  inertial  force 

(S.Cif{  ) ‘/^ry-)  is  the  external  force 
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Then  *;  Structural  force 


- (m.) 


^ 0/  "f  • 

But  TLOif,  = and  X is  obtained  in  terms  of  . 

Therefore  ^ eliminated  from  equation  (111;.  mt) 


6.  It  is  clear  that  the  generalized  acceleration  factor  does 
not  appear  as  such. 

7.  Using  the  "double"  trapezoid  deriyed  above,  the  right  side 

of  equation  (110)  is  already  scaled  relative  to  Pmay*  In  the 
computation  of  Table  3,  the  ordinate  of  the  trapezoid  will  be 
multiplied  by  Pmair  so  that  the  structural  force 

1/2  Airplane  Weight 

is  finally  non->dlmensionallzed  in  units  of  1/2  airplane  weight. 

Step  V.  Vibratory  Acceleration  Response: 

Table  3 computes  the  vibratory  response  by  the  desk-calculator 
method. 

Step  VI.  Total  Acceleration  Tima  History: 

The  final  dimensionless  result  of  ^ 

One-half  weight  of  airplane 

is  ccmiposed  and  oompar^d  with  experimental  results  in  Figure  ll. 

The  exqperimental  force  normal  to  the  strut  at  the  axle  was  computed 
from  the  measured  drag  link  load  by  taking  a summation  of  mcmtents 
about  the  top  of  the  strut  assrxning  the  oleo  to  be  fUUy  extended. 


* Of  course  the  differeritial  equation  (5)  could  be  solved  this 
time  for  X rather  than  x , but  it  is  convenient  to  utilize  the 
methods  of  Section  II. 
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ITlALLr.COMpUTABLE:  CONSTANT3  IZ.  In  AB.  e"  Sin  Cut 
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